import sympy as sym

x,y,z=sym.symbols('x,y,z')
f1=2*x**2+3*x+1
print(f1.subs(x,2))
print(sym.factor(f1))
sym.plot(f1)

f2=f1*(x+1)
print(f2)
print(sym.expand(f2))

result=sym.solve(f1)
print(result)
